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fS| . Abstract 
> 

Q\ • Products of random matrices associated to one-dimensional random media satisfy 

a central limit theorem assuring convergence to a gaussian centered at the Lyapunov 
exponent. The hypothesis of single parameter scaling states that its variance is equal 
to the Lyapunov exponent. We settle discussions about its validity for a wide class of 
models by proving that, away from anomalies, single parameter scaling holds to lowest 
order perturbation theory in the disorder strength. However, it is generically violated at 
higher order. This is explicitely exhibited for the Anderson model. 
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1 Introduction and main result 



One-dimensional quantum systems with a single channel can very efficiently be described by 
2x2 transfer matrices. In a disordered medium, the transfer matrices are chosen to be random. 
The one-dimensional Anderson model is the proto-type for this class of models. The most 
important physical phenomena in these random media is localization due to multiple coherent 
wave scattering. It goes along with positivity of the Lyapunov exponents associated with 
products of the random matrices and the Lyapunov exponent is then interpreted as the inverse 
localization length of the system. Moreover, the fluctuations around this asymptotic behavior 
are gaussian. More precisely, Anderson, Thouless, Abraham and Fisher [ATAFJ stated that 
the Landauer conductance follows asymptotically (in the system size) a log-normal distribution 
centered at the Lyapunov exponent. As pointed out by Johnston and Kunz |JKj . this was 
a rediscovery of a mathematical result by Tutubalin |Tutj (refined by Le Page |LePj ) . The 
paradigm of single parameter scaling |AALR[ lATAFj is then that there is only one parameter 
describing the asymptotic behavior of the random system. For a one- dimensional model this 
means that the Lyapunov exponent and the variance of the gaussian should be in some relation 
and in fact simply be equal (ATAFj . The validity of single parameter scaling in this sense has 
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been analyzed in various particular situations ATAFj ISAJ| ICRS| IDLA| ISTj . The main result 



of the present work can roughly be resumed as follows: in a wide class of one-dimensional 
random models single parameter scaling is valid only perturbatively in a weak disorder regime 
and never holds in a strict sense (exceptional parameter values excluded). 

For this purpose, we study a general class of one-parameter families of random transfer 
matrices exhibiting a so-called critical energy. This parameter is the effective size of the ran- 
domness. In the Anderson model, it is the coupling constant of the disordered potential while 
in the random dimer model DWP it is the distance in energy from the critical energy. We then 
develop a rigorous perturbation theory in this parameter. In case of the Lyapunov exponent, 
we do not appeal to the random phase approximation ATAFl ISA Jj . but rather show that phase 
correlations give a contribution to lowest order perturbation theory (which, however, vanishes 
for the Anderson model). This generalizes arguments of |Tho| IPF[ IJSSj . On the other hand, 
the perturbation theory for the variance is new to our best knowledge. One has to sum up 
the phase correlation decay by using adequate counter terms. This rigorous analysis is made 
possible by a result of Le Page (see Section EJ) . Furthermore, we extend the techniques of [SS 
in order to calculate the scaling exponent of the expectation value of the Landauer conductance 
(this is sometimes also called a generalized Lyapunov exponent). This generalizes results of 
iMoll. 



Comparing the coefficients to lowest (namely second) order perturbation theory away from 
Kappus-Wegner type anomalies |KWj , we obtain that the Lyapunov exponent and the variance 
are equal while the scaling exponent of the averaged Landauer conductance is twice this value, 
just as predicted by |ATAFj . Calculation of higher orders is possible, but cumbersome in 
general. For the Anderson model it becomes feasible and is carried out in Section El We 
obtain that the next (namely forth) order contributions are not the same for the Lyapunov 
exponent and the variance and as a consequence they are not equal (but close) in the regime of 
weak disorder. For the regime of strong disorder, even large discrepancies have been observed 
numerically |SAJj . Deviations from single parameter scaling even to lowest order were exhibited 
at the band center of the Anderson model, which is the prime example of a Kappus-Wegner 
anomaly |STj . The Lloyd model analyzed in |DLA| does not fit in our framework because there 
the random variables do not have finite moments. The example of the Anderson model also 
allows to show that there does not exist a universal analytic function expressing the variance 
in terms of the Lyapunov exponent. 

After this brief introduction, let us describe our results more precisely. The transfer matrices 
are supposed to be elements of the following subgroup of the general linear group Gl(2, C): 

U(l, 1) = {T6 Mat 2x2 (C) | T*JT = J} , J = [l ~q ) • (*) 
Using the conjugation 

c* jc = i r , c = -i I ! \ \ , r ( 1 



V2 V 1 - 1 J ' V - 1 

one sees that U(l, 1) is also isomorphic to the subgroup of matrices T G Gl(2, C) satisfying 
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T*TT = T. This representation appears in some applications, but for our purposes it is more 
convenient to work with (JI} because it contains the standard rotation matrix and the real 
subgroup SL(2,R). 

We will study families (T A -) AgK/Te2 of transfer matrices in U(l, 1) depending on a ran- 
dom variable a in some probability space (E, p) as well as a real coupling parameter A. The 
dependence on A is supposed to be smooth. 

Definition 1 The value A = is a critical point of the family (Tx tU ) xm,a£T: if for all a, <r' 6 S: 

(i) [T 0>a ,T 0!a r] = , (ii) |Tr(T 0i(T )| < 2 . (2) 

Critical points appear in many applications like in the Anderson model and the random 
dimer model |DWP( I.TSS[ ISedj , but also continuous random Schrodinger operators where the 
transfer matrix is calculated from a single-site S'-matrix |KSj . Condition (i) assures that there is 
no non-commutativity at A = (even though the matrices may be random), while by condition 
(ii) the matrices T 0)Cr are conjugated to rotations so that there is no a priori hyperbolicity in 
the system. The example of the Anderson model is studied in more detail in Section El 

Associated to a given semi-infinite code {a n ) n >i is a sequence of matrices (T\,a n )n>i- Codes 
are random and chosen independently according to the product law p® n . Averaging w.r.t. p® N 
will be denoted by E. We will suppose that all up to the 5th moment of p exist. In order to 
shorten notations, we will also write T\^ n for T\^ n . Of interest is the asymptotic behavior of 
the random products 

JV 
n=l 

It is first of all characterized by the Lyapunov exponent 

7 (A) = lim i-Elog(||T A (iV)||) . (3) 

N^oo iv 

The central limit theorem for products of random matrices now states that 

1 'log(||T A (iV) e ||)-iV7(A)) GWa) (4) 



where G^x) is the centered Gaussian law of variance cx(A) and the convergence is in distribution 
independently of the initial unit vector e. This was first proven by Tutubalin under the hy- 
pothesis that the measure p has a density |Tutj . Le Page then proved it for arbitrary measures 
P |LePj . Both proofs can be found in |BLj . As already discussed above, the single parameter 
scaling assumption is the equality er(A) = 7(A) [ATAF, CRS . Apart from the Lyapunov expo- 
nent 7(A) and the variance <r(A), we are going to analyse the growth exponent of the average 
of the Landauer conductance defined by 

7(A) = lim -^log(ETr(T A (A0*T A (A0)) • (5) 
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It follows immediately from Jensen's inequality that 7(A) > 7(A). Our main results are resumed 
in the following: 

Theorem 1 Introduce the phase r\ a G [0, 27r) by cos(r] a ) = iTr(T 0jCT ). Suppose E(e 2yr ' tT ) 7^ 1 
for j = 1,2,3. Then there is a constant D > 0, given in equation (}2*2|) below, such that near 
the critical point, 

7(A) = DX 2 + C(A 3 ) , 7(A) = 2D X 2 + C(A 3 ) . (6) 
// moreover D > 0, then 

a(X) = DX 2 + C(A 3 ) . (7) 

In Section |H we also give criteria insuring that D > 0. The remainder of the paper contains 
the proof of this theorem as well as an analysis of higher orders for the Anderson model. 

Acknowledgments: This work profited from financial support of the SFB 288. H. S.-B. would 
like to thank T. Kottos for pointing out references jDLAj and [ST . 



2 Normal form of transfer matrices near critical point 

If T G U(l, 1), then T* JT = J implies that det(T) = e 2 ^ for some f G [0, vr). Hence e~*T G 
SU(1, 1) = {T G U(l, 1)| det(T) = 1}. This means U(l, 1) = U(l) x SU(1, 1). Note also that 
Tr(T) = det(T)Tr(T) for T G U(l, 1) so that Tr(T) G R for T G SU(1, 1). One easily verifies 
that this implies SU(1, 1) = SL(2,R) is a real subgroup of G1(2,C). Its Lie algebra is well 
known: 



sl(2, R) 



a b 
c —a 



a, b, c G 



hi 



The eigenvalues of T G SL(2,R) always come in pairs k,1/k where k = Tr(T)/2 + 
lyjl — Tr(T) 2 /4 and are hence complex conjugate of each other if | Tr(T) | < 2 (which is the case 
of all transfer matrices near a critical energy). In this situation, the matrix diag(/t, 1/k) is not 

(cosiT/l — sin (n) \ 
. ) ( / \ with k 
sm(?7) cos(r]) J 

is in SL(2,R) and can be attained with an adequate conjugation. Hence for any T G U(l, 1) 
with \Tr(T)\ < 2, there exists M G SL(2,R) and rj such that MTM~ l = e^R v . 

Let us now consider the family (Tx^AeiveE i n U(l, 1) satisfying (J2J). Because they commute 
and have a trace less than 2 for A = 0, they can simultaneously be conjugated to a rotation at 
that point. Using a Taylor expansion in A, we therefore obtain the following: 

MTx^M- 1 = e^ (A) R Va exp (XP a + X 2 Q a + 0(X 3 )) . (8) 
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Here e 2 ^^ = det(T A>(7 ) so that \P a + X 2 Q a E sl(2,R) for all A G R. In particular, Tr(P a ) = 
Tr (Qcr) = 0. The constant C in © only depends on the rotation angles r\ a and on P a through 
the constant 



0a 



(v\P v \v) 



(9) 



v = — ( 1 

Note that = (v\P*\v), \p a \ 2 = \Tr(\P a \ 2 + P 2 ) and that P 2 = (P CT *) 2 is a multiple of the 
identity. Moreover, v and v are the eigenvectors of all rotations R a . 

Given a n , let us now denote the associated random phases, rotations and perturbations by 
£ n (A), Tj n , R n , P„, Qui thus suppressing the dependence on the random variable a n . 

In order to use the normal form (jHJ) for the calculation of the Lyapunov and Landauer 
exponents, let us insert M~ X M in between each pair of transfer matrices. As it only gives 
boundary contributions, one may also insert an M to the left and an M~ l to the right. Hence, 



7 (A) = lim ^Elog 

N— >oo iv 



A' 



II MT \nM~ l 



n=l 



as well as a similar formula for 7(A). Now it is clear that there is no use in carrying along 
the phases e %£,a ^ in (jSJ) if one is interested in calculating 7(A) and 7(A). Therefore we may set 
from now on £ CT (A) = 0. This is equivalent to supposing that T\ >a £ SL(2, R). Moreover, one 
may factor out the subgroup {1, —1}, namely even work with the projection of MT\^ a M~ x into 
PSL(2, R) = SL(2, R)/{1, -1}. 



3 Further preliminaries and notations 

Unit vectors in M? will be denoted by: 

Each transfer matrix T\ a induces an action on unit vectors via 

MT Ka M- l e e 
e s xA e) - \\ MTxnM -i ee \\ ■ 

Using the vector v of (JOJ), this is equivalent to 



MTx^M^ee 



(v 




ee) 


(v 


MTx^M- 1 


eg) 



(10) 



(12) 



Now given an initial condition 9 , a random sequence of phases 9 n associated to a code (cr n ) ngN 
is iteratively defined by 



9 n — S\,a n {9n-l) 



(13) 
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A probability measure v on S 1 is called invariant for this random dynamical system if 



dv{9)f{9) 



dv(9)E a f(S x<a (9)) , feCiS 1 ) 



where E CT denotes the average w.r.t. p. Due to a theorem of Furstenberg [BLJ v exists and is 
unique whenever the Lyapunov exponent is positive. Positivity of the Lyapunov exponent for 
A 7^ is guaranteed by condition (iii) of Definition Q (we leave it to the reader to verify that 
the subgroup generated by the transfer matrices T Xa is then non-compact so that Furstenberg's 
criterium is satisfied |BL| ). In the sequel, E v will mean averaging w.r.t. v as well as the whole 
code (cr n ) neN . 

Next let us turn to the Lyapunov exponent. According to |BI4 A. III. 3. 4] it is given by 



7(A) = lim — E log 



N 



\\ MT x>n M~\ 



n=l 



(14) 



As one can also insert the invariant measure, it is also given by the so-called Furstenberg 
formula: 



7 (A) = J dv{9)E a log {\\MT Xi<J M- l e e \\) . 
Finally, one can use the random phase dynamics (|T3~j) in order to rewrite (fTf)) as 



7(A) 



1 - 

N— >oo iv ' 



n=l 



(15) 



4 Asymptotics of the Lyapunov exponent 



Let us introduce the random variable 



7 n = log (|| MTx, n M~ 1 ee n _ 1 1|) . 
Then the results cited in the previous section imply 

N 



7 (A) 



lim , 



In 



lim E (7„ 



71=1 



Our first aim is to derive a perturbative formula for 7„. Replacing 
orthogonal and expanding the logarithm shows 



(16) 



(17) 



using that R n is 



In = ^KJP n |eO+y(e^^ (18) 
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where we have used that (P 2 )* = P 2 is a multiple of the identity and set P n = P n + P* and 
Qn = Qn + Qn- Now for an Y T e Mat 2X 20R), one has 

(e e \T\e e ) = ~Tr(T) + e 2lfl ) . 

Hence using the definition © and the remark following it, we deduce 



In = y|/5„| 2 + ^e^ n e 2 ^--y fie 4 *- 1 + y (U|(|P n | 2 + Qn)H e 3 *- 1 ) + 0(A 3 ) . 

(19) 

The so-called random phase approximation consists in supposing that the angles Q n -\ are 
distributed according the Lebesgue measure (i.e. v is the Lebesgue measure). Then only the 
non-oscillatory term in (|T9*j) would contribute so that one would get 7(A) = |AE cr (|/5 -| 2 ) + (9(A 3 ). 
In general, however, this is erroneous. Replacing (fT9*|) into (fTKj) . one has to calculate the 
following oscillatory sums (as in JSSJ). 



Lemma 1 For j = 1,2, set 

N-l 



2jiO n 

N 



n=0 

Suppose E(j (e 2jM?CT ) ^ 1 for j = 1,2. Then 

h(N) = 1 _ g £ g ( ^ ) J + 0(\\ iV-) , UN) = 0(X, iV-) . 
Proof. It follows from (JT2D and flSJ) that 

((l + AP^Kl + AP^) +U[X} - (20) 

In particular, this implies 

so that, replacing this in each term, one obtains 

Ij(N) = E CT (e 2 ^) Ij(N) + O^N- 1 ) . 

The induction hypothesis therefore implies that Ij(N) = 0(X, N' 1 ). In order to calculate the 
contribution of 0(A) to h(N), let us expand (J2Uj) . Some algebra shows that 



e 2i6»„ 



; 2 ^" +e - l) - Ae**- 1 ' (e 2j9 "- 1 /3 n -2(t;|P n |t;) -e^" 1 /^) + £>(A 2 ) . (21) 



From the last three terms, those containing still an oscillatory factor e 2l0n or e il9n will not 
contribute to leading order A due to the above. Thus we deduce 
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h(N) = E a (e 2 ^) h(N) + AE CT (e 2%Va j3 a ) + 0(X 2 , N' 1 ) . 

This implies the result. □ 

As lirriAT^oc Ij(N) = E u (e 2lje ), the lemma shows that E u (e 210 ) = 0(A) and calculates the 
lowest order contribution. If (3 a is centered (as for the Anderson model), one even has E u (e 2ie ) = 
0{\ 2 ). This shows in particular how far the invariant measure u is away from the Lebesgue 
measure and that phase correlations are indeed present. 

Now we can proceed with the calculation of 7(A). Carrying out the algebra shows 7(A) = 
D X 2 + C(A 3 ) where 

c ^ W) + fe ( g^pp ). (22) 

The second summand is due to phase correlations. It is important, e.g., in the random polymer 
model [JSS . In the Anderson model treated in Section El phase correlations only contribute to 
the forth order in A. Because 7(A) > 0, the coefficient D defined by (J22j) has to be non-negative. 
More precisely, we prove: 

Proposition 1 Suppose E CT (e 2 *^) 7^ 1. Then D is always non-negative. D vanishes if and 
only if one of the following two mutually excluding cases occurs: 

(i) Both e 2w?CT and (3 a are p-a.s. constant. 

(ii) E CT (e 2M?CT ) = and f3 a is a constant multiple of 1 — e 2%r,a . 

Proof. Assume first that e 2vn " is a.s. constant. Then E fT (e 2lJ?<T ) = e 2jr?<T 7^ 1 and it follows 
from - e^)' 1 = 1/2 that 2 D = E CT (|/3 CT | 2 ) - \E a ((3 a )\ 2 . By the Cauchy-Schwarz inequality, 
D > and D = if and only if (3 a is a.s. constant. Now let us assume that e 2%r>a is not a.s. 
constant so that |E CT (e 2M?CT ) | < 1. The proposition then follows from the following lemma by 
setting H = L 2 (p), ^1 = 1 and V> 2 = e 2w?CT . □ 

Lemma 2 Let ipi and ip2 be two linearly independent unit vectors in a Hilbert space TC, implying 
K^i 1^2) I < 1- Then the quadratic form 

on positive semi- definite. It is positive definite if and only if (ip\\ip2) 7^ 0. If {ip\\fp2) — 0; 
i/iera Q(ip) = i/ and onfo/ if ip is a multiple of ipi — 

Proof. Let K, be the two dimensional subspace of TC spanned by ipi and ip2 and /C -1 its 
orthogonal complement. For an arbitrary vector ip, we write ip — ip' + ip" with ij}' G K and 
•0" G JC ± . Then Q(V>) = Q(V ;/ ) + HV'"!! 2 - Hence it suffices to assume ip — ip' e JC, i.e. we may 
restrict Q to /C. If we introduce the orthonormal basis 

ei = ipi , e 2 = == === (ip 2 ~ A ) 

V 1 - \mm)\ 2 
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in /C, then with respect to this orthonormal basis the quadratic form Q on /C is represented by 
the 2x2 matrix 

/ l+K^lv-a)! 2 yi-|(^i|fe)| 2 \ 

(l-W>l|</>2»(l-W>2|V>l» 1-W>2|V>1> 



Q 



1 



/ 



\ l-(Vl|^2> 

Now its trace is obviously strictly positive, while its determinant satisfies 

2|(^i|^)| 2 



detQ 



> o . 



Hence both eigenvalues of Q are strictly positive whenever ip\ and ip2 are not orthogonal, while 
one eigenvalue is strictly positive and the other one zero if they are orthogonal. For this zero 

eigenvalue the corresponding eigenvector of Q is proportional to the vector 



concludes the proof of the lemma. 



This 

□ 



5 Perturbation theory for the variance 

Using telescoping as in (JHj), the variance is given by 

cr(A) = lim — E ^ (7n - 7) (7k - 7) 



n,k=l 
N 



lim 



N-n 



^E E ^ 7^-7 2 + 2^( 7 n7 

n=l m=l 



n+m 7 



where we wrote 7 for 7(A) for notational simplicity. Let us show that the sum over m is 
convergent even if N — > 00. To this aim, we denote by E n the expectation over all a m with 
m > n such that with the previous notation Ei = E. From the remark following Proposition 121 
in Section El it follows that E n ( 7n+m ) converges exponentially fast to 7 for m — > 00. Moreover, 
the summands of the sum over n converge in expectation so that, if v is the unique invariant 
measure (integration w.r.t. 9q), the variance is given by 



(7(A) 



7l 2 - 7 2 + 2 (7i7m - 7 2 ) 

m=2 



E„ 



7l 2 - 7 2 + 27 1 E 2 (7m -7) 

m=2 



• (23) 



Hence we need to evaluate the sums appearing in the following lemma. Its proof is deferred to 
Section El 



Lemma 3 Suppose E CT (e 2ljl)a ) ^ 1 for j = 1,2,3 and D > 0. Th 



en 



E 2 (7m - 7) 



A 3?e 



m=2 



E CT (/3 CT ) 
1 - E a (e 21 ^ 



+ 0(X 2 
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As 7 = 0(X 2 ), the first summand 7 2 gives no contribution. Therefore we obtain 



o-(A) = E„ 



2 2 A E -(/3 CT ) 2 <e +m) 

71 + Ke 1-E ff (e**) 71 



+ 0(A 3 ) . 



Let us calculate the first contribution supposing that v is the invariant measure. Then extract- 
ing the linear coefficient in A from (fTTj|) and using Lemma Q 



E ^ 2 > = & f e "£ 7 » 

n=l 

= ^ lim -i E„ £ (A^+^e- 2 ^) 2 + 0(A 3 



Similarly 



This implies 



4 TV^oo iV 

n=l 

^ E CT (|/? CT | 2 ) + 0(A 3 ) . 



E,( 7l e 2 ^)) = ^E CT (^e 2 ^) + 0(\ 2 ) 



<t(A) = 7 (A) + 0(\ 3 ) . 



6 Estimates on the correlation decay 

The main purpose of this section is to prove Lemma El We choose a slightly more general 
formulation, however, allowing to treat also other quantities. The first aim is to show how the 
estimates of |BL| Theorem A.V.2.5] (due to Le Page |LePp can be made quantitative when 
combined with Lemma ^ Throughout this section, we suppose that D > and, for sake of 
notational simplicity, A > 0. Let us introduce a distance on S 1 by 



8(9, ip) = yl-(e e |e^,) 2 = ||e e Ae^||, 

where the second norm is in A 2 IR 2 . For a G (0,1], the space of Holder continuous func- 
tions Cai^S 1 ) is given by the continuous functions / G C(S 1 ) with finite Holder norm 
max{||/|| 00 ,m a (/)}, where 

m 1/(0) - /Ml 

m a (f) = sup — — . 

e^esi 0(9, ip) a 

Then we define S, V : C a (S\) — > C a (Si) by 

(Sf)(9) = E(f(S x , a (9))) , (Vf)(9) = J 'duty) fW . 
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Here V should be thought of as the projection on the constant function with value given by 
the average w.r.t. the unique invariant measure. Both S and V depend on A. 

Proposition 2 There exists a constant c, d > such that for a = d X, 

\\(S N -V)(f)\\ a < ||/|Ue- cA3Ar . 

Proof. Let us introduce ( : SL(2,R) x S 1 x S 1 -> R by 

where iSr : S 1 — ► S 1 is defined as in (fTT|) by es T (o) = Teg/\\Teg\\. Then £ is a cocycle, namely 
it satisfies 

C(T'T, (9,ip)) = C(T',(S T (e),Sr(m + t(T,{0,i/>)) . 

Moreover, one has 

C(T,(fi,^)) = logf ||T ini^ii^i ll ) = -log(IM)-log(||Te^||), 
V||Te || IITe^H \\e e A e^\\J 

because ||A 2 Tee A e e+ ^\\ = 1. Therefore ^(MT^M -1 , (6 1 , ip))\ < c A and the cocycle property 
implies 

\({T x {n),(9,i>))\ < Cl \n. 
Furthermore, invoking Lemma Q shows that 

E(C(T A (n), (M))) < -2DA 2 n+c 2 (A 3 n + A) , 

where D > is the coefficient given in (|2"2"j). Defining the sequence of angles ^ n as in ([TSjl. but 
with initial condition ipo = ip, one can now infer from e x < 1 + x + x 2 & x \/2 that 

E ( WW) - 1 + «E(C(T A (n),(e,V))) + ^ 2 E(C(T A (n),(e^))Vl^^))l) 

< 1 - 2 L> a A 2 n + c 2 a (A 3 n + A) + - c\ (aAn) 2 e Cl oAn . 

2 

Now one has on the one hand, 

and, furthermore, using the invariance of z/ and then 5(0, ip) < 1, 
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\\(S n -V)(f)l 



sup 

ees 1 



dv{^) E(f{d n )-fU> n )) 



< m a (f) sup / dvfy) E 



Therefore, we can deduce 



\\(S n -V)(f)\\ a < m a (f) (l-2aD\ 2 n + c 2 a(X 3 n + X) + cl(aXn) 2 e CiaXn ) . 

Now we choose n — [c 3 /A] + 1 (as usual [x] denotes the integer part of x G R) and a = c 4 A 
with adequate c 3 , c 4 > and note that m a i(f) < m a (f) for a' < a, so that for some c 5 > 

\\(S n -V)(f)\\ a < \\f\\ a (l-c 5 X 2 ), n = 

Finally, we note that S N — V = (S n — V) N ' n because SV = VS = V so that iterating the last 
inequality completes the proof. □ 

When applied to the function f{6) = E cr (log(||MTA iCr M~ 1 e6i||)) (which is Holder continuous 
for any a) this estimate directly implies that E(7„) converges exponentially fast to 7 when 
n — > 00. Furthermore, using m a (f) < \\dgfWoo, one now gets: 

Corollary 1 There is a constant c > such that 

00 

E (/(0tn)-E„(/(0))) < cmax{||/|| 0O ,||^/|| oo } A" 3 . 

m=l 




We will now consider an algebra A of functions which are analytic on some neighborhood 
of {0} x S 1 C C 2 and have a power series of the form 

fc>0 \l\<k k>0 

with complex coefficients F^i. The only elementary fact needed here is that F-G G A whenever 
F, G G A. We will mainly be interested in the values on R x 5* 1 and there also write F(X, 6) = 
F(X, 2 2td ). Let us give two examples for elements in A. One is the dynamics defined in (fTTj) . 
namely (A, 6) h-> e 2%< ^ Sx ^^~ 6 ^ is a function in A because it is a quotient of analytic functions in 
(A, 9) of the form (J24j) . Taking its jth power, one obtains the representation 

e 2 V s xA e) = e 2^ + „ CT ) [l + ^A^c^e 2 ^] , (25) 

y k>i \i\<k J 

for adequate complex coefficients c^(cr). Comparing with (jUJ), we see that c\ x {a) = —f3 a , 
c\ (a) = 2 (v\P a \v) and ci_ 1 (cr) = -/3 CT . 
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Our second example is the function F(X, 9) = E cr (log(||MTA j(J M 1 ee> 1 1 ) ) - With this function, 
one has E m (7 m ) = F(X,9 m -i) and 7(A) = E„(F(A, #)). This function appears in LemmaE! 

Generalizing Lemma El we are led to evaluate (perturbatively in A) the summed up corre- 
lation decay for an arbitrary function F G A: 

00 

Cor(F)(A) = E 2 ^2(F(X,9 m )-E v (F(X,9))) . 

m=l 

From Corollary ^ follows the a priori estimate |Cor(F)(A)| < C A~ 3 . For simplicity, let us 
now calculate Cor(F)(A) to order A, namely discard terms of order 0(X 2 ). This is the sit- 
uation covered by Lemma H Consider G(X,9) = F(X,9) - ELo Afe = so that 
\\deG{\, . )|U = C(A 5 ). Combined with Corollary d it follows: 



Cor(F)(A) = E 2 J2 G(A,0 m ) + J]A fc F fc (^ m )-J]A fc E 1 ,(F fe (0))-E i ,(G(A,0)) 

m=l \ k=0 k=0 

4 00 

= ^A fe E 2 ^2(F k (9 m )-E v (F k (9))) + 0(X 2 ) 



k=0 m=l 

4 00 



k=l \l\<k m=l 

The appearing sum over m is finite (again due to Corollary |TJ) and, moreover, we can calculate 
its value perturbatively. 



Lemma 4 Let 



Jj = E 2 J2 i e2tjdm ~ Vu{e 2tjd )) , J E Z 



m=l 



Suppose E„(e 2 «'-) / 1 for j = 1, . . . ,4. TTien J 4 = 0(\- 2 ), J s = 0(A" 1 ), h = O(l) and 

- 71 = 1-E.W) e2 ' < ' ,0+ '" , + ° (A) ' (26) 

Proof. We use (I25j) in order to express e 2 *- 76 '" 1 in terms of e 2y m_1 , but truncate the expansion 
at C(A' K "). Using again Corollary [TJ we deduce 
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J, 



m=l 



K 



1 + E Afe E E ^>)) e 



2iW„ 



k=l 



\<k 



K 



2%w 



0(X 



K-2\ 



- Ecr (e 2 ^) E„ e 2 ^ 1 + £ A fc £ E ^>)) e2 
\ \ fc=i |i|<fe 

E CT (e 2tjv ' y ) Jj + e 2 ^ 9o+ ^ + E^e 2 ^) ^ A fe ^ E ff (c^(cx)) J j+ , + 0(A^ 2 , A) 

fc=i \i\<k 

K 



E CT (e 2 ^' 



,2ii(£/ +7?l) 



+ E. (e^j^A^E.^^))^ 



fc=i 



|<jfc 



0(A K - 2 ,A) , 



where C(A^" 2 , A) = 0(A^ 2 ) + 0(A) and in the second equality we used the fact E u (e 2ljd ) = 
(9(A) (due to Lemma Q). As we know that Jo = and that the a priori estimate Jj = 0(X~ 3 ) 
holds, this calculation shows that Jj = C(A" 2 ) if E CT (e 2lj ^) ^ 0. The estimate J 3 = cjA" 1 ) 
now follows by choosing K = 1 and replacing J 4 = 0(A~ 2 ) and J 2 = 0(X~ 2 ). The same way 
one deduces J2 = 0(X^ X ) and J\ = 0{X^ 1 ). Choosing K = 2, a similarly argument shows 
in the next step J 2 = 0(1) and Ji = 0(1). In order to establish (j2EI), let us choose K — 3. 
The first term in the last square bracket gives the desired contribution, while the sum is 0(A) 
because J 2 = 0(1), J 3 = O^A" 1 ) and J 3 = 0(A" 2 ). □ 

In order to calculate Cor(F)(A), let us recall that J = and J_j = Jj. Thus 

Cor(F)(A) = XF ljl J 1 + AFi_i Ji + 0(A 2 ) , 

with Ji given by LemmaH For the function F(X, 9) = E a (log(||MT Aj(T M- 1 e e ||)), equation (flU 
implies F^i = ^E a (/3 a ) = F1-1. This gives Lemma 01 



7 Some identities linked to the adjoint representation 

The adjoint representation of SL(2, K.) on its Lie algebra is defined by Adr(t) = TtT~ l , t G 
sl(2,R). It leaves invariant the quadratic form q(t,s) = ^Tr(ts) of signature (2, 1) (note that 
q(t, t) = - det(t)). A basis B = {61, b 2 , b 3 } of sl(2, R) is given by 

61 = ( -1 ) ' 62 = ( 1 I ) ' 63 = ( 1 'o 1 ) ' 

It is orthonormal w.r.t. the scalarproduct (s\t) = |Tr(s*t). Moreover, denoting the coordinate 
map w.r.t. this basis by K B : sl(2, M) — > M. 3 and the standard scalar product in M 3 also by 
(x I y) for x, y G M 3 , we have 
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q(t,s) = (K B {t)\T 2)1 \K B {s)) , r 2jl = 

Finally, let us set Ad? = K B ; Ad T (i^ B ) _1 . This means (Ad|) i|fc = (6 i |Ad T (6 fc )). Note that 
J = b 3 G SL(2,R) and that Adf = -r 2 ,i. Thus TVT = J implies (Adf )*r 2>1 Ad T = r 2>1 . 
Hence Ad B is an element of the Lorentz group S0(2, 1) given by all A G Mat3 X 3(M) satis- 
fying v4*i\iv4 = r 2j i. As Ad B = Ad B x = 1, the adjoint representation gives an isomorphism 
PSL(2, R) = SO(2, 1). Important for the sequel is that the eigenvalues fix, /i 2 , /X3 of A G SO(2, 1) 
satisfy \i\\i2^z = ±1. One of these eigenvalues, say /X3, must be real, while the other two may 
either be real as well or be complex conjugates of each other. 

Next it follows from a short calcuation that for T = f « ] G SL(2,R) 

\ c a J 

(ad + bc db — dc dc + bd 

dc-ao \{d 2 -b 2 -c 2 + d 2 ) \{d 2 -b 2 + c 2 -a 2 ) 
dc + dh \{d 2 + h 2 - c 2 - d 2 ) \(d 2 + l 2 + c 2 + d 2 ) 

In particular, for a rotation R v by rj we get 

cos(2r]) — sin(2r/) 
Ad^ = [ sin(2??) cos(2r?) 
1 



The eigenvalues are e 2M? , e 2lv , 1 with respective eigenvectors V\ = (ei — ze 2 )/v2, t? 2 = (e*i + 
ze 2 )/\/2, w 3 = e 3 . 

We will also need the adjoint representation of the normal form (jHJ: 

Adp exp(XP+A 2 Q) = Ad R ^l + Aadp + A 2 ad Q + y (ad P ) 2 + £>(A 3 )^ , (27) 

where ad f (s) = [i, s] for t,s G sl(2, R). Let us write out more explicit formulas in the represen- 
tation w.r.t. the basis 8. For P = ( a b J G sl(2,K), 



c —a 

6-c o + c 

ad P = [ c-b -2a | , (2i 

b + c -2a 



and 



46c -2a(b + c) 2a(c - 6) 

(adp) 2 = ( -2a(6 + c) 4a 2 - (c - 6) 2 c 2 - 6 2 

2a(6 - c) b 2 - c 2 (b + c) 2 + 4a 2 

Finally let V : C 2 <8> C 2 — > C 2 <E> C 2 be the permutation operator 
readily be checked that 
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-l 



Multiplying this identity from the left by T <g> J and from the right by 1 Cg> T J, one gets for 
T G SL(2,R), 

T T'P = \ ( - J® J + ± detfe) (Ad ?U , W JJ,) . (30) 

This is useful for the calculation of the Landauer conductance because 

\\Tw\\ 2 = (w\Tri(T®T t 'P)\w) , 

where Tri is the partial trace over the first component of C 2 (g) C 2 . Replacing (J3U|) . the first 
term vanishes because Tr(J) = 0. Moreover, Tr(6fcJ) = —2 5^3 so that 

3 

||Tw|| 2 = ^(Ad£) 3 j(-l)det(&j) (w\Jbj\w) . 
Let us define a map g : C 2 — > M 3 by 



WlW 2 + ^1^2 

12 L„ 12 „,. _ / Wl 



( det(6i) 




J6i|w) \ 




det(6 2 ) 




Jb 2 \w) 


-( 


V det(6 3 ) 


(w 


Jb 3 \w) J 





//if) = - 1 act^o) \«V»2|<'7 | = i F2I - Fil I ' w ir 



\W 1 + f 2 



Then one has 



||TH| 2 = (e 3 \A4\g{w)) . (31) 
Since g(t>) = e 3 , this implies in particular 

||T W || 2 = <e 3 |Ad*|e 3 ) . (32) 

8 Calculation of the averaged Landauer resistance 

The Landauer resistance p\(N) of a system of length N is defined by 

p x (N) = E (Tr(|T A (iV)| 2 )) = 2E (\\T x (N)v\\ 2 ) , 

with v as in (|9|) . Because \T\(N)\ 2 is positive and has unit determinant, one has p\(N) > 2 for 
all N. Using the identity (JH2J) and then the representation property Adf T = Adf Ad^ iteratively, 
the expectation value appearing in the Landauer resistance can readily be calculated: 

Px (N) = 2(e,\(E(Adl)) N \e 3 ) , 
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so that 7(A) = liniAr^oc log(p\(N)) / {2N) . Replacing the normal form (JHJ) (recall that the phases 
disappear right away in the definition of the Landauer conductance): 

Px (N) = 2<(Ad«- 1 ) t e 3 | (EAd B RacMXPa+x2Q ^ +om) ) N \Ad B M e 3 ) . 

Next we need to do (non-degenerate, in A) perturbation theory of the eigenvalues pi, p 2 , P3 
of EAd| CT exp(AP CT +A 2 Q CT +o(A3)) ; which according to (J27J) is up to C(A 3 ) given by 

E(Adl) + AE(Adlad^) + A 2 E(Adladgj + yE(Adl(ad|j 2 ) , 

Let us note that the eigenvalues of E(Ad^ ) are pi = E(e 2ir?CT ), p 2 = E(e~ 2ir>,T ), p 3 — 1 
with the eigenvectors v\, v 2 and v 3 respectively. Therefore, unless r\ a is independent of a, the 
matrix E(Ad^ ) is not an element of SO(2, 1) and two of its (complex conjugate) eigenvalues 
are strictly within the unit circle. 

We first focus on the eigenvalue p 3 . Because (e3|adp|e3) = 0, one has p 3 = 1 + 0(X 2 ) and 
the eigenvector is t7 3 + 0(X 2 ). Second order perturbation theory now shows 



A 2 



H 3 = 1 + - (v,\E(Ad^(adf,J 2 )\v 3 ) 
+ A 2 (^|E(Adlad^)- 



E(Ad 



;i-|^>(^|)E(Adiad?J|i%> + 0(A 3 



Recalling (b + c) 2 + 4a 2 = 4 \(3 a \ 2 and using (v 3 \&dpjvi) = -^(b + c + 2ta) = i\f2(3 a , we deduce 



/'3 



1 + 2A 2 



E(|/3 C 



+ 2 5Re 



1l1] a 



+ (D{\ 3 



(33) 



1 - E(e 2 ^) 

Now let us analyse the eigenvalues p\ and \xi- If | E(e 2tr,a )\ < 1 (as for the dimer model), they 
are strictly within the unit disc and remain there also for A sufficiently small. As p\(N) > 2, 
we conclude that p 3 > 1 (this also implies D > 0, like in Proposition [TJ because otherwise 
Pa (AO — > for — > 00). It follows that the scaling exponent of the Landauer resistance 
defined in (JHJ) is given solely by p 3 , namely 7(A) = \{p 3 — 1) + C*(A 3 ) so that, when comparing 
with (J22I), 7(A) = 27(A) + C(A 3 ) as claimed in Theorem □ 

In the case where 7] a = rj independently of a (as in the Anderson model), one has D > 0. 
First order perturbation theory (in A) shows that pi^ move along the the unit circle, while in 
second order they lie inside of the unit circle. Thus the same argument as above applies to 
deduce Theorem [TJ 



9 Higher orders for the Anderson model 

This section serves two purposes: it provides an example to which the general theory applies 
and we moreover outline (algebra is left to the reader) how it can be extended to calculate the 
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next higher order in perturbation theory. The implications of the results have already been 
discussed in the introduction. 

The one- dimensional Anderson model is a random Jacobi matrix given by the finite differ- 
ence equation 

-VVfl + Xv n 1p n -tpn-i = Elp n . 

Here \E\ < 2 is a fixed energy and v n are centered real i.i.d. random variables with finite 
moments. This equation is rewritten as usual using transfer matrices: 

1pn+l \ T ( Ipn \ rp f XVn-E -1 

As above, we also write v a for one of the random variables such that v n = v an and T A (Tn = T A n . 
For the basis change to the normal form of the transfer matrix T\ a , let us introduce 



E = -2cos(fc) , M 
It is then a matter of computation to verify 



sin(fc) 
^/sin(fc) V -cos(k) 1 



MT^M- 1 = R k (l + AP CT ) , P a = - -^tt ( ? n 

sm(/cj \ 1 U 

Comparing with (jHJ, we see that the rotation i?^ by the angle k is not random in this example, 
that P a is nilpotent so that exp(AP CT ) = 1 + XP a and that Q a = 0. Furthermore 

^ = > = ■ 

One has (3 2 a = -\(3 a \ 2 , ft = \(3 a \ A and 

(e e \P a \e e ) = p a (e™ - e~ 2ie ) , (e e \ \P a \ 2 \e e ) = \(3 a \ 2 (2 + e™ + e" 2 * 9 ) . (34) 

As E(/3 CT ) = 0, one can immediately deduce from (|2*2*|) the well-known formula |Tho[ IPF1 ILucj 
for the Lyapunov exponent, namely 7(A) = §A 2 E(\(3 a \ 2 ) + £>(A 3 ) = A 2 + 0(X 3 ) as long 

as e 2y ' fc 7^ 1 for j = 1,2 (the latter condition excludes band edges and the Kappus-Wegner 
anomaly at the band center E = 0). 

In order to calculate the 4th order of the Lyapunov exponent and the variance, we need 
higher order expansions than ()21j) and ()18|1 . We will assume e 2lJ 7^ 1 for j = 1, ... ,4. After 
some algebra, 



{v\R k (l + XP a )\e e ) 
(v\R k (l + XP a )\e e ) 



(35) 



O 2i(o+k) 



Xf3 a {e™ + 2 + e~ 2ie ) + A 2 /3 2 (e 4 ^ + 3e 2 ^ + 3 



+ C(A 3 
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and, using 

lx , a (9) = log(||(l + AP CT 
= 



A ',^ e - \ 2 f3l ( h 4ld + + i J + A 3 /? 3 Qe 6ie + e 4i6 + e 2ie 



(36) 



+ a 4 - e™ - - e™ - ^ + G(A 5 ) 



Using an argument similar to Lemma ^ we deduce from (|35|) and its square that 

A 2 E(|/3 CT | 2 



E„(e 



2i8\ 



+ C(A d 



El/ ( e 4 ^ = X * E{lPa D + 0(A 3 ) 



^ _ e -2ik 1 v / ' "V / j_ _ e -4ifc 

Moreover, E^e 64 ") = C(A 2 ) and E^(e 8 * e ) = C(A 2 ). Following the argument of Section EJ we 
therefore obtain from (|36|) and 3?e(l — e^) -1 = | 

7(A) = | A 2 E(|^| 2 ) + A 4 (^E(|/? CT | 2 ) 2 - jE(|/3 CT | 4 )) + G(A 5 ) . 



This coincides with the fourth order contribution obtained in |Lucj by using complex energy 
Dyson-Schmidt variables. 

In order to calculate the variance according to formula (|23|). one first needs to go through 
the arguments of Section One finds 



E 2 ( e2l9m ~ E,(e 2i9 )) 



O 2i8 



-2ik 



m=l 



- (l-A/?i(e 2i9o +2 + e- 2 ^)) + C(A 2 



a similar expression for the correlation sum of e , and that E 2 £^=2(7?" ~~ t) ^ s u P t° ^(A 
equal to 



A 2 E(|/3 ff | 2 ) 9£e 



e 2 ^(l - A/3i(e 2je ° + 2 + e" 2 ^ )) 1 e 4ie °(l - 2\f3 1 (e 2te ° + 2 + e~ 2%e °)) 

+ 2 



g-2«fc _ ]_ 

Finally, after having squared (|3*5jl. 

E,( 7l 2 )-7 2 = ^A 2 E(|/3 ff | 2 ) 

and 



| A 4 E(|/3 CT | 4 ) - iA 4 E(|/3 CT | 2 ) 2 + 0(A 5 ) 



Ej2 7l E 2 £)(7m-7) = ^A 4 E(|^| 2 ) 2 + 0(X 5 ) . 



m=2 



Combining these results according to (J23j) . we obtain 
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*(A) = | A 2 E(|/? CT | 2 ) + A 4 Qe(|&| 2 ) 2 - |E(|/3 CT | 4 )^) + 0(\ 5 ) . 

Therefore we see that Lyapunov exponent and variance are only equal to lowest order in per- 
turbation theory. Finally let us also argue that there cannot exist a universal analytic function 
/ such that a = f("f). Indeed, if f(x) = fix + f 2 x 2 + 0(x 3 ), then cr(A) = /(7(A)) for the 
Anderson model implies in order A 2 that f± — 1, but in order A 4 there is already a problem due 
to the prefactors of E(|/3 CT | 4 ). 
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